Mixed moment wave function for magnetic heavy fermion compounds 
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We propose a variational wave function for the ground state of the magnetic heavy fermion 
(HF) systems, in which both the Kondo and the RKKY interactions are variationally incorporated 
and the local f-orbital state exists as a linear combination of a full local moment state and a fully 
compensated state (mixed moment state). We describe the mechanism for the mixed moment 
ground state based on the large-N treatment of the Kondo lattice Hamiltonian added with RKKY 
interaction. With the mixed moment ground state we can explain several puzzling experiments in 
magnetic HF compounds such as a small value of local moment, coexistence of the antiferromagnetic 
(AFM) and the paramagnetic (PM) phases, local quantum criticality, etc. 



I. INTRODUCTION 



As a conventional wisdom, one often employs the 
Kondo mechanism to understand the HF behavior, in 
which the conduction electrons make bound states with 
the local f-orbital moments in a coherent fashion to result 
in fermionic quasiparticles integrating the local f-orbital 
moments with it. While it is still unclear how the single 
site Kondo mechanism can be generalized to the peri- 
odic lattice system, at least the above line of thinking 
seems to provide the correct energy scale, i.e. the coher- 
ent energy scale of forming heavy renormalized fermionic 
quasiparticles 1 . However, this phenomcnological picture 
doesn't explain many experimental observations of HF 
systems, such as, strong temperature dependence of spe- 
cific heat C(T) and bulk susceptibility x(T) below the 
coherent temperature T co h, the coexistence of the mag- 
netism (often antiferromagnetism (AFM)) and the heavy 
quasiparticle in some compounds, and superconductiv- 
ity in some other compounds. Those compounds which 
displays both HF properties as well as the magnetism 
(either magnetic fluctuations or the magnetic long range 
order(LRO)) are conveniently called "magnetic HF" by 
experimentalists. Then again the common wisdom for 
this class of HF systems is the competition between the 
Kondo coupling between the local f-orbitals and the con- 
duction electrons and the RKKY interaction among the 
local f-orbitals themselves 2 . And at qualitative level the 
result of this competition is that the ground state should 
be either pure HF state (all local moment is completely 
compensated by conduction electrons) or pure magnetic 
state (no compensation of the local moment) depending 
on which energy scale is larger. However the experimen- 
tal situation is that these magnetic HF compounds not 
only show the coexistence of the local moment and the 
heavy quasiparticle but also almost always the size of 
the local moment is much smaller than its ionic limit 
value. Therefore, it is clear that a minimum requirement 
to understand experimental observations is that the true 
ground state wave function for the magnetic HF com- 
pounds should simultaneously possess both the heavy 



quasiparticle and local moment characters. 

In this paper, we propose a variational wave function 
for the ground state of the magnetic HF systems, in which 
the local f-orbital wave function exists in the mixed states 
(a linear combination of the full local moment state and 
the fully compensated state), and we call it "mixed mo- 
ment state" after the "mixed valence state" 3 . In the fol- 
lowing sections, we show how this mixed moment state 
can be realized in the Kondo lattice model using the 
large-N scheme. This part is rather qualitative based 
on account of the energetics of the competition between 
different couplings. Once this mixed moment state is phe- 
nomenologically accepted, then this ground state prop- 
erties can immediately provide explanations for the most 
puzzling questions of magnetic HF compounds: small un- 
compensated local moment, the coexistence of the AFM 
and PM domains (the PM domain can be even a su- 
perconducting phase), and as the most interesting pos- 
sibility, the local quantum criticality coinciding with the 
magnetic criticality and its relation with the Fermi sur- 
face (FS) fluctuations, etc. 



II. THE MIXED MOMENT STATE 

We start with the Kondo lattice Hamiltonian (KLH). 

k,rn i,m 
, J Kondo 



N 



E 3 - . Irkky sr^ q c 



where cf k m ,Ck, m is the conduction electron creation and 
annihilation operators and f\i m ,fi, m arc the local f- 
orbital operators, and the spin operators are the large 
N generalization of SU(2) spin Si — f\i W rf i m > and 
Si = cf i m fc i m > for the f-orbital and conduction elec- 
tron spins, respectively. The above Hamiltonian is the 
usual Kondo lattice Hamiltonian extended with the large 
N spin degeneracy (m = 1, N) and the f-orbital occu- 
pancy is constrained with n{ = £ /ti,m/»,m = ^<7o 4 
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representing the large on-site Coulomb interaction be- 
tween the localized f-orbitals at the same site. Finally the 
RKKY interaction term between the f-orbitals is added 5 . 
The above Hamiltonian has been studied by numerous 
authors using various techniques 6 . Without the RKKY 
term, the various approximation (mainly the approxima- 
tion treating the constraint n{ — Nq ) provides a solu- 
tion of heavy renormalized coherent band(s) to explain 
the HF phenomena. When the RKKY term is added as in 
Eq.(l) not only any approximate solution becomes more 
complicated but more importantly the question of the 
correctness of treating the constraint becomes crucial be- 
cause the fate of the competition/interplay between the 
Kondo and the RKKY interactions is largely unknown 
to determine the true ground state. There are a few 
studies of the two impurity version of this model . All 
these studies indicate that there is a critical ratio of in- 
teractions {J Kondo I Irk ky) from which the ground state 
flows away toward either a magnetically correlated state 
or the Kondo singlet state. While we learn from these 
studies that there is an interesting competition between 
the Kondo and RKKY interactions, a naive extension of 
this picture to the lattice system is not warranted. 

As discussed in the introduction, motivated by experi- 
mental observations, we construct a new variational wave 
function for the ground state of the above Hamiltonian. 
We start with the f-orbital wave function as follows. 



\fi,m >= \/l - Oi\fi >m >iti +Va\fi,m >loc 



(2) 



The above expression is designed to indicate that the 
f-orbital state in the ground state can be a superposi- 
tion of two qualitatively different states; one is the itin- 
erant state (|/i, m >iti) which makes a coherent singlet 
state with the conduction electrons via renormalized hy- 
bridization through Kondo coupling 6 and the other is the 
full local moment state (|/j, m >i oc ) which remains intact 
from forming Kondo singlet or heavy quasiparticle but 
couples only as a local magnetic moment with other elec- 
trons. With the above Ansatz for the f-orbital state, we 
diagonalize the Hamiltonian Eq.(l) in large N approxi- 
mation. The ground state should determine the value of 
a to minimized the total energy of the Hamiltonian. This 
wave function, if it is realized as a true ground state, dis- 
plays the mixed moment state in a basically same manner 
as the mixed valence state 3 , and the origin of these mixed 
states is the local Coulomb interaction. 

Assuming the variational Ansatz of Eq.(2), the Kondo 
coupling term in Eq.(l) is decomposed into two parts 
depending on which component of f-orbital states is in- 
volved. 

J Kondo \ ~* ci ^ /i \ J Kondo 

Si-ai = (l-a)- 



N 



N 



J2 /ti,mCi,mCti, m ' Urn' + a ^j^ ^ Si ■ ff { (3) 



In the above equation, we rewrite the Kondo term cou- 
pled with the itinerant component of the f-orbital by the 



four fcrmionic expression often employed in the large-N 
treatment of Kondo coupling 6 . The second term repre- 
sents the local moment part of the f-orbital state, which 
remains as a local moment and does not participate in 
forming the heavy quasiparticles. Now the Hamiltonian 
can be written as, 

H = H Kondo + Hmag, (4) 

HKondo — J]] ek,mC\k,m c k,m + e / J]] /ti,m/i,m (5) 
k,m i : m 

+ (l-a) JK ™ do J2 fU, m Ci,mtfi,m'fi,m'' 
i,m,m' 

rr J Kondo ST^ Z -. . ilRKKY \^ X Z / g n 

H mag = a — 2^Si-ai+a N2 Si ■ Sj (6) 



HKondo part is well studied by many authors and in 
particular the large N technique is convenient to produce 
the heavy renormalized quasiparticle band(s). Switch- 
ing on the Kondo coupling {J Kondo), the ground state 
of HKondo lowers its energy compared to J Kondo = 
state (AE Ko ndo)- H mag part also should gain energy by 
developing a magnetic correlation, but the first term of 
H mag , by coupling with the conduction electrons (c elec- 
trons), also increases the kinetic energy of the HKondo- 
First, we consider HKondo- The only difference from the 
previous studies is that the Kondo coupling JKondo is 
replaced by (1 — a)JKondo and also the f-orbital occu- 
pancy constraint is modified as n{ — (1 — a)Nqo (ac- 
cordingly the Fermi surface volume changes with a). 
Introducing the Stratonovich-Hubbard decoupling field 

<f>0 = (l- aJ^Em < fU,m C i,m >, H K ondo IS diagO- 

nalized and the renormalized e/ is fixed to satisfy the con- 
straint n{ = (1 — a)Nqo. Following Read et al. (Ref.[2]), 
the Kondo energy gain AEKondo is written as, 



AE K ondo = N P v4>l H^f/D] - - a ) + 



N 



Jk(1-o) 



(7) 



where D is the typical conduction band width and p a 
is the density of states (DOS) of it. The self-consistent 
equations are as follows. 



Npo4>l- = (l-a) 
ef = D exp [ 



PqJk{^ - a) ' 



(8) 
(9) 



Using the above equations, we find the Kondo energy 
gain per site is the following. 



AE Ko ndo = -(1 - a)Dexp[ 



-1 



(10) 



The system gains the above energy by compensating a 
fraction (1 — a) of the local moments via Kondo screening 
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and the conduction band becomes a renormalized heavy 
band(s). Now let us consider H mag . The RKKY term 
would gain the following energy per site with the maxi- 
mum polarization < Si >= S z = ^r-- 

AE RKKY = ~a 2l -^Sl ~ 0(N°). (11) 

The first term of H mag (Eq.(6)) now forces the AFM 
coupling of the conduction electrons (c electrons) with 
a staggered field hf or a staggered energy level Aj = 
a^-S z . With this staggered energy level for c-electrons, 
the renormalized heavy band(s) should further develop 
spin density wave (SDW) ordering and it will cost the 
kinetic energy increase as follows. 



AE kin = \p HF ul F ^ + ^0z 



Q) 2 



(12) 



where phf is the DOS of the renormalized band at chem- 
ical potential and Uk is the Bogoliubov coefficient of Ck 
component of the renormalized band operators which di- 
agonalize Eq.(5). v F is the Fermi velocity and Q is the 
SDW ordering vector 8 . Finally the total energy gain of 
H mag with both the local magnetic and SDW orderings 
is as follows. 



AE. 



mag 



'Phf Wfc^A 2 



2 IRKKY c2 



(13) 



Now collecting all the energy gain and loss, the total 
energy difference by the Kondo and RKKY couplings is 
written as 



AE tota l = AE 



Kondo ~t~ AE^in "i" AE ma g 
-1 



= —(1 — a)Z?exp[ 



PoJKondo(.^ — a) 



(14) 



(15) 



where A^ = -\p HF ^(^^ 
AW 



ll W LS l and 

v%{ P k H /_ Q) A ^rSi. Note that A^ < 0, 

A^ > and since S z ~ O(N) they are all 0(N°). Ellipsis 
indicates the higher order terms (0(a 6 ) etc.) 

Treating A^ and A^ as phenomenological parame- 
ters, in Fig.(l) we show the schematic total energy gain 
AE to tai(oi) of Eq.(15) for the three representative cases. 
There is always a local minimum with a finite value of 
a (due to A^ < 0, A^ > ). However, when this lo- 
cal minimum is not a global minimum as in Fig.l.a, the 
ground state is a pure HF state (a = 0, no magnetic HF) 
with fully screened local f-moments. In contrast, when 
this local minimum with a finite a becomes a global min- 
imum as depicted in Fig.l.b., then the ground state is 
a "mixed moment" state. If a fine tuning occurs such 
as using pressure, magnetic fields, chemical substitution, 
etc., then two minima can be degenerate with a ^ and 
a = as in Fig. I.e. In this case, as the tuning parameter 
changes, the ground state goes through a first order tran- 
sition accompanying with a jump of magnetization and 



FS volume. Finally, the above discussions are mean field 
level and therefore a a smaller than a critical value would 
not develop a true long range order (LRO) of local mo- 
ments due to fluctuations; it means that the nonmagnetic 
phase near the magnetic phase boundary should have a 
small but finite a value. 



III. NEUTRON DIFFRACTION (ND) AND NMR 
EXPERIMENTS IN URU 2 SI 2 



In our mixed moment state, the ground state is a su- 
perposition of the unscreened full local moments and the 
Kondo renormalized itinerant band(s). We first consider 
the consequence of this ground state for the measure- 
ment of the local moment by neutron diffraction (ND) 
experiment. At any given time, a fraction (a) of U sites 
have a full local moment and a fraction (1 — a) of U 
sites are completely quenched by the conduction elec- 
trons. Now the unquenched local moment is magneti- 
cally screened by the opposite SDW moment developed 
by the renormalized itinerant band. This incompletely 
screened magnetic moment is the effective local moment 
(S e J f = S z -m z , where m z = n co „ d (|) - n cond {i)). Now 
the ND should see the ensemble averaged size of the lo- 
cal moment, i.e. p, e ff = gfisS z ^ x a. This value can 
be tiny (p, e ff ~ 0.03^b for URu 2 Si2) while S z ^ and m z 
are not necessarily tiny. 

A recent NMR experiment by K. Matsuda et al. 10 re- 
ports another puzzling data on magnetic moment. In this 
experiment, 29 Si NMR spectra clearly shows that 29 Si 
sites see both the AFM ordered phase and the PM phase 
below the magnetic transition temperature To and these 
authors interpret the data with a spatial inhomogeneous 
mixture of AFM and PM domains. From the estimated 
volume fraction of AFM domain, it is concluded that the 
actual size of local moment is at least 0.3^.b/U, an order 
of magnitude larger than the estimation from the neutron 
Bragg peak intensity 11 . In our picture, the coexistence of 
AFM and PM is not of a static inhomogeneous domain 
structure, but at any given time the 29 Si nuclei should 
see a fraction (a) of and a fraction (1 — a) of zero 
moment from U sites. From this experiment we can read 
the size of SDW moment by S e J^ = S z —m z = 0.3/ g and 
S z = p-para/^B = 1-2/g. Combining ND and NMR data, 
Me// ~ 0.03^s = (0.3fiB x a), a ~ 0.1 is estimated for 
URu 2 Si 2 for ambient pressure. Also the observation of a 
distribution of local effective fields at 29 Si sites from 29 Si 
NMR by O.O. Bernal et al. 12 is a natural consequence of 
the local moment fluctuations of nearby U sites between 
zero and a finite value. 

A strong pressure dependence of the local moment 
size from ND by H. Amitsuka et al. 13 (from 0.017^b to 
0.25^b) and much slower increase of the magnetic order- 
ing temperature (T m ) is not inconsistent with our model: 
the measured local moment size is (U full local moment 
— SDW moment) x« and in zeroth approximation T m 
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is ~ a 2 (Eq.(6)) assuming the RKKY coupling Irkky 
constant 9 . Applying pressure should change a (increas- 
ing in URu2Si2) and T m increases, but the increase of the 
effective local moment size should be a more complicate 
function of a. 



IV. FERMI SURFACE VOLUME AND C(T) JUMP 

In the mixed moment state, the total carrier density 
in the renormalized itinerant band is given as n to t — 
n c + Nqo(l — a) (n c is the original conduction electron 
density). Therefore, the FS volume changes by how much 
of the f-electrons participates to form the coherent heavy 
quasiparticle band(s) through the Kondo coupling. As 
a result the fluctuations of the local moment weight a 
is intimately coupled to the fluctuations of FS and any 
phase transition of one part should trigger a transition 
of the other part. One consequence is an enhanced spe- 
cific heat jump AC(T) at the second order transition as 
commonly observed in many magnetic HF compounds. 
In general, the specific heat jump at the second order 
transition temperature can be expressed as follows. 



AC mag /T 



d 2 F 

ACsc/T ■ 



a i 



bv 



dT 
OA 2 
dT 



a-2 



dT 



dT 



(16) 
(17) 



The first case (m; oc = aS z ) can be applied to URu2Si2, 
which is well known for its extremely tiny local moment 
(fj, e ff ~ 0.03/is) from ND but displaying a huge jump 
(AC/T ~ 300m J / K 2 mol) at the magnetic transition 
T = 17.5K. To obtain some quantitative scales, first note 
that T = 17. 5K indicates the RKKY coupling should 
be large enough to justify the relatively high transition 
temperature with a small averaged local moment value 
(aS z ), and it also means that J Kondo is not a small value 9 
and indeed the experimental Tr- ~ 60K . Thus the aver- 
age internal fields induced by the local moment ordering 
is given by h\ = _g_ J Ka^da g z (^he first term in Eq.(6)), 
which is not so tiny value despite the acclaimed small 
effective local moment. Then the itinerant renormalized 
band develops a SDW ordering over a part of FS simul- 
taneously at T — To and it can provide the major contri- 
bution to the specific heat jump (in the mean field theory 
«2 ~ Phf)- The second case (Eq.(17)) is for the SC tran- 
sition (for example, CeCoIn 5 14 ). The deformation of FS 
due to the SC should trigger a concomitant change in a. 



V. COINCIDENCE OF LOCAL QUANTUM 
CRITIC ALITY AND THE MAGNETIC 
CRITIC ALITY 

Many magnetic HF compounds exhibit non-Fermi liq- 
uid (NFL) behaviors in resistivity, specific heat, neutron 
scattering, susceptibility etc 16 . Often these behaviors 



coincide with the magnetic quantum transition (mostly 
AFM transition). However, theoretically any quantum 
critical fluctuations, which is spatially correlated, is not 
sufficient to explain the NFL behaviors. Therefore, what 
is naturally required is a "local quantum criticality" not 
only for NFL HF compounds but perhaps also for the 
high temperature superconducting materials. 

There are a few theoretical proposals 17 ' 18 of local quan- 
tum criticality in the strongly correlated electron system 
with different ideas. Our "mixed moment" ground state 
can provide a natural mechanism for the local quantum 
criticality, in particular, for NFL HF compounds. The 
first term of Eq.(6) provides a coupling between the con- 
duction electrons (ex,) and local spins (Si), and the lo- 
cal spins develop its own dynamics through the RKKY 
interaction, which can go through the magnetic quan- 
tum criticality by tuning external parameters. This is a 
two component spin-fermion model recently proposed to 
study the superconductivity of CcMIns 15 . Now in addi- 
tion to the spatially correlated magnetic fluctuations, if 
the moment weight fluctuations < SctiSctj > is allowed 
in higher order corrections, this fluctuations absorbs any 
value of momentum exchange just like impurities as far 
as the typical energy of this fluctuations is much higher 
than the typical low energy scale of the spin fluctuations 
< SiSj >. Then all conduction electrons around FS can 
be scattered off each other by a critical magnetic fluctu- 
ations with a help of the local moment weight fluctua- 
tions. This process is depicted diagramatically in Fig. 2. 
Therefore, in this picture the spatially correlated mag- 
netic quantum criticality becomes a local quantum criti- 
cality with a help of the local moment weight fluctuations. 



VI. CONCLUSION 

To summarize, we propose a new variational wave func- 
tion for the magnetic HF which represents a mixed mo- 
ment ground state. We use a mean field analysis in large 
N limit of the KLH added with RKKY interaction and 
show that this is indeed a generic ground state for a cer- 
tain ratio of the couplings of Kondo and RKKY. Then 
we show that this mixed moment ground state can im- 
mediately provide natural explanations for the most puz- 
zling observations in magnetic HF compounds ( URu2Si2, 
UPt3, CeMIns, CeCus-zAu^, etc.) such as a tiny or- 
dered magnetic moment, large specific heat jump, a coex- 
istence of AFM and PM phases, a distribution of internal 
fields, and most interestingly the local quantum critical- 
ity coinciding with the spatial magnetic criticality. More 
quantitative calculations for each of these properties will 
be reported in separate publication. 
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Science and Engineering Foundation (KOSEF) through 
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FIG. 1. Energy gains as a function of a for the repre- 
sentative cases of different ratios of Kondo and RKKY cou- 
plings: AEtot(a) (solid line), AEKondo(a) (thin solid line), 
AE mag (a) + AE kin (a) (dotted line). 




FIG. 2. Feynman diagram for the vertex of electron elec- 
tron interaction mediated by magnetic fluctuations. 
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